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Abstract 

We consider the discrete Poincare constant, which relates the variance of a function 
to the expected square of its finite difference. We give an expHcit bound on the 
Poincare constant of ultra log-concave random variables in terms of their first two 
moments, and discuss how this bound relates to calculations performed by other 
authors. 
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1 Introduction 

In this article, we consider the finite difference operator A defined by {Ag){x) = g{x + 
l)—g{x) for functions (yf : N = {0, 1, . . .} — M. We discuss probability mass functions Px 
supported on N for which Ag being small in L'^{Px) implies that g is close to constant 
(in the same sense). We make the following definition: 

Definition 1.1 Define the Poincare (spectral gap) constant of a random variable X 
with mass function Px supported on N to be 

R - sun ^^-oPx{x)g{xr 
rCx — sup 
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where the supremum is taken over g in the set 

{oo "\ f oo ^ oo 

J2Px{x)g{x) = n $^Px(x)/(x) < oo fl lj2Px{x){Ag){xr > 
x=0 ) I x=Q ) L x=0 

The motivation for this article is to show that the ultra log-concave (ULC) random 
variables (see Definition 12.31 below) defined by Pemantle [12] and Liggett [11] have finite 
Poincare constant. In fact we can give a bound in terms of the first two moments. 

Theorem 1.2 

(i) If X is an ultra log-concave random variable of degree oo then it has Poincare 
constant Rx satisfying 



Rx <EX + ^ + ^^ + EX - Vai X. 

(a) If X is an ultra log-concave random variable of degree n then it has Poincare 
constant Rx satisfying 



1 EX 1 EX 
fl.<EA- + 5- — + ^j + E.Y-Var.Y-— . 

Proof See Section [2] below. □ 

The plan of the paper is as follows. In Section O we give some techniques for calculating 
Poincare constants. In particular, Lemmas 12.11 and 12.21 give conditions on the tails of 
the random variable X under which the Poincare constant Rx is bounded, allowing us 
to prove Theorem 11.21 

The idea of proving such discrete Poincare inequalities is not a new one, with previous 
authors to consider the problem including Bobkov and Gotze pj, Bobkov and Ledoux 
[2], CacouUos and Papathanasiou [3], Chen and Lou |6], Klaasen [H], Prakasa Rao and 
Sreehari [13] and Srivastava and Sreehari 1141 . In Section [3] we discuss how our results 



compare with theirs. 

In Section H] we correct a mistake in Part (ix) of Theorem 2 of [3] . 



2 Calculating Poincare constants 

We require a direct equivalent of Theorem 1 of Borovkov and Utev [3] . 
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Lemma 2.1 Let X have probability mass function supported on N of the form 

Px{x) = aPi{x) + (1 - a)P2{x), 
where < a < 1. Suppose that for some xq and c> 0, the inequalities 



^ (x - Xo)Px{x) < cPi{y) for y > xq, (1) 

x=y+l 

y 

^(x - Xq)Px{,x) < cPiiy) for y < Xq. (2) 



x=0 

hold. Then 



Rx < c/a. 



Proof Equations ([H) and ([2]) are based on those given by Borovkov and Utev [3], and 
are applied to the proof given by Klaasen [9]. For the sake of completeness, we give 
the whole argument here. The key is to consider Klaasen's kernel function, given by 
Equation (2.17) of [9]: 

y) = I( [xoj <y <x)-I{x<y < [xqJ ) - (xq - [a^oj Mv = [xo\). (3) 

For any given integer x, by considering the cases {x < xq}, {x > xq} and {x = xq} 
separately, we deduce that for any function h: 

°° [ - T.l=t^ Hy) + ( l^ol - Xo)h{ [xo\ ) for X < Xq, 

^X{x,y)h{y) = I ([xoj - xo)/i(La;oJ) forx=[xoj, (4) 

y=o ( Y.l~^^ h{y) + ([xoj - Xo)h{[xo\) for x > Xq. 



Taking h = (Ag) we deduce that J2'^=oX{x,y){Ag){y) = g{x)-g*, where g* = g{[xo\) + 
(A5f)(LxoJ)(xo — [xo\). In particular, taking h{y) = 1 we deduce that ^^o^(-^'^) ~ 
(x — xo). Observe that by Cauchy-Schwarz: 

{g{x) - g*f = (f^xix, y) {Ag) {y)\ < (x - Xo) x{x, y) {Ag) (y) A . (5) 
\y=o J \y=0 ) 

Further, by assumption. Equation ([3]) implies that for any y, 

oo 

Y^x{^,y)Px{x){x - Xq) <cPx{y). (6) 

a;=0 

[Equation ([6]) follows by considering three cases separately, using Equations ([I]) and ([2]): 



3 



(i) For y > \_xq\ the left-hand side becomes Yl^=y+i Px{x){x — Xq) < cPiiy) by ([T]). 

(ii) For y = [xqJ the left-hand side becomes (1 — xq + [xo\) Yl'^=y+i Px{x){x — xq) < 
(1 - xo + [xo\)cPi{y) < cPi{y) by ([I]), since for any xq, < 1 - xq + [xqJ < 1. 

(iii) For y < [xqJ the left-hand side becomes Yl^x=a-Px{.x){x — Xq) < cPi{y) by 
This means that (with the reversal of order of summation justified by Fubini): 

oo oo 

Y,Px{x)g{xr < Y.Px{x){g{x)-cf 

oo / oo \ 

< 5^Px(x)(x-Xo) 5^x(a;,l/)(A<7)(l/)M (7) 
x=o \y=o ) 



oo 



< 



Y,{^9){yf\^Px{x){x-x^)x{x,y) 

y=0 \x=0 



< cY,{^9){yfPM (8) 

3/=0 
oo 

-Y,{l^g){yfPx{y). 



< 



and the result holds, since for all y, the Pi{y) < Pi{y) + (1/a — l)P2{y) = Px{y)lot- 
Note that Equation ([71) follows by Equation ([5]), and Equation (IHl) follows by Equation 
(0). □ 

Note that Cacoullos and Papathanasiou [1] prove a similar result, but under the more 
restrictive condition that Xq = E^. 

Lemma 2.2 

(i) If there exists C such that 

xPx{x) . { >C forx<C, 

then the Poincare constant Rx < C . 

(ii) If for some n there exists D such that 



An),. ._ ^Px{x) / >D/{l-D) forx<Dn, 

{n-x + l)Px{x-l) \<D/{l-D) forx>Dn, ^'^^ 

then the Poincare constant Rx < Dn 
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Proof (i) If Equation © holds then taking xq = C, for y > xq: 



Px{x){x-xo)<C J2 Pxix-1)-C J2 Px{x)=CPx{y), 

x=y+l x=y+l x=y+l 

and for y < Xq: 

y y y 

J2Px{x){xo - x) < cY^Pxix) - cY^Pxix - 1) = CPxiy), 

x=0 x=0 x=0 

SO that Equations ([I]) and ([2]) hold with ci = C, so by Lemma [2.11 Rx < C. 



If Equation (fTOj) holds, then for x > Dn then we can rearrange in a way suggested 
by Klaasen's work to obtain 

Px{x){x - Dn) < D{Px{x -l){n- x + 1) - Px{x){n - x)), 

so taking xq = Dn, for y > xq, since the sum collapses: 

oo oo 

Px{x){x - Dn) = D Y Px{x-l){n~x + l)-Px{x){n-x) = D{n-y)Px{y), 

x=y+l x=y+l 

which is < DnPx{y)- A similar argument works for y < xq with the signs reversed, so 
that Equations ([T]) and ^ hold with ci = Dn, so by Lemma I^TTl Rx < Dn. □ 

Lemmas 12.11 and 12.21 imply that the Poincare constant is finite for the following classes: 

Definition 2.3 A random variable X is ultra log-concave of degree n, or ULC(n) for 

short, if X is supported on {0, 1, . . . ,n} and the function Px{x)/(2) is log-concave, or 
equivalently 

(n). . xPxjx) . . . . 

p\- \x) = -, , ^ , IS non-mcreasmq m x. 

^ ' {n-x^\)Px{x-\) 

A random variable X is ultra log-concave of degree oo, or ULC(cxd) for short, if the 
function Px{x)x\ is log-concave, or equivalently 

xPx{x) . . , , 

Pxix) = ; r IS non-mcreasmq in x. 

Pxix-1) 

These classes were introduced by Pemantle [12j and Liggett flT] in order to capture 
properties of negative dependence. Liggett remarks that if a random variable is ULC(n) 
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then it is also ULC(n + 1), and proves in Theorem 2 of [Ll\ that if X and Y are 
independent and ULC(n) and ULC(m) respectively then their convolution is ULC(?t,+ 
m). Hence the sum of n independent Bernoulli random variables has the ULC(?2) 
property. Clearly the Poisson random variables are ULC(oo). 

Proof of Theorem 11.21 If X is ULC(oo) then px{x) = xPx{x) /Px{x — 1) is non- 
increasing in X, and r(x) = x is strictly increasing in x, so Equation must hold for 
some C < oo, so by Lemma 1^:^ Rx < C. We can analyse Equation ^ further to 
deduce bounds on C. Considering the cases {x < C} and {x > C} separately, we know 
that for each x the summand is negative, so deduce that the sum 

> y - 1) f T^T^^ -C\{x-C)= - C(2EX + 1) + C\ 

Treating this as a quadratic in C, we know that this means that C G [C_, C+], for some 
roots C±. Hence, we need simply calculate the value of C+. 

If X is ULC(?2) then p'^\x) = xPx{x) / {Px{x — l)(n — x + 1)) is non-increasing in x, 
and r(x) = x/{n — x) is strictly increasing for x < n, so Equation ( ITOi) must hold for 



some D < oo, so by Lemma 12. ^(ii) , Rx < Dn. We can analyse Equation (ITO!) in a 
similar way to deduce bounds on D. Considering the cases {x < Dn} and {x > Dn} 
separately, we deduce that the sum 

> Ed - mPA^ - - ^ + A'it- X + 1) - T^) - 

= EX^ - Dn{2'EX + 1 - EX/n) + D'^n^. 
Treating this as a quadratic in C = Dn, wc know that 



^ 1 EX /I EX EX 2 EX 2 ^ ^ 

C<EX + + W- + EX-VarX ^ + ^ 11 

2 2n V 4 2n n 4n^ ' ^ ' 

and the sum of the last two terms is negative, so they can be removed to give a neater 
bound. □ 



Remark 2.4 The upper bounds given in Theorem \l.S\ are indeed real-valued, since the 
variances of ultra log-concave variables are constrained in terms of their means A. Specif- 
ically, for X in ULC(r;,), since we can view it as the covariance of a decreasing and 
increasing function (with respect to probability mass function {n — x)Px{x) /{n — X)) the 
sum 

{n-x)Px{x) r Pxjx + l){x + 1) X_ 

n — X \ Px{x){n — x) n — X^ 



E 



X 
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is negative, by Chebyshev's rearrangement lemma. Rearranging, this tells us thatE,X{X— 
1) ^ X'^{n — l)/n, or that EX — Var {X) > }? jn. (The corresponding result is well known 
forn = 00). Substituting in the bounds of Theorem \1.2\ we deduce that the argument of 
the square roots are positive, and indeed the same is true for Equation ( flij) . 



3 Comparison with other results 

We now discuss how our results compare with known bounds. First we mention a simple 
lower bound, the direct equivalent of Theorem 2(vi) of [3J: 

Lemma 3.1 For any random variable X with finite variance: Rx > Var {X). 

Proof Consider the function g{x) = {x — EX), so that Ag{x) = 1. We know that 
E^(X)2 = VarX and E{Ag){Xf = 1, so that Rx = supg{Eg{Xf)/{E{Ag){Xf) > 
Var (X)/l. □ 

We briefly discuss how our results differ from those of other authors to have considered 
discrete Poincare constants. 

Bobkov and Gotze [Ij give necessary and sufficient conditions for the discrete Poincare 
constant to be finite. However, as we shall see, their methods can significantly overesti- 
mate the exact value of the constant. 

Theorem 3.2 ([Ij) Given random variable X with probability mass function Px , define 

where Fx{x) = YTy=QPx{y)- Bobkov and Gotze' s work implies that 

C{Px) 



C{Px) <Rx< 



Px{0) ■ 



Proof The left-hand side follows by considering gx{y) = ^{y < x) — Fx{x) for each x. 
Then g^iyf = I{y < x){l - 2Fx{x)) + Fx{x)\ so Egl{Y) = Fx{x){l - Fx{x)), and 
{Ag,){y) = I{y = x) so E{Ag,){Yy = Px{x). 

The proof of the right-hand side is more complicated, but the result is stated as the 
third displayed equation on P.274 of [T]. □ 
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Our Lemma 12.11 was designed to generalise the calculation of the Poincare constant of a 
Poisson random variable, given in Table 2.2 of Klaasen's paper [U] and Proposition 3.5 
of Cacoullos and Papathanasiou [5] in the same year. Lemma 12.11 reduces to give the 
values they calculated in this case. 

Example 3.3 Writing Il\{x) = e~^X^/x\ for the Poisson mass function, we know that 
xll\{x) /Ilx{x — 1) = A for all x. Hence Equation ^ holds with C = \. Lemma \2.^ij\ 
implies that Ru^ < A. Combining this with Lemma l3J\ we deduce that 

Rn, = A. (13) 

Taking x = in Equation H^) we obtain that Cijlx) > 1 — 11^(0) = 1 — exp(— A), so 
Bobkov and Gotze's Theorem \3.2\ implies an upper bound that is at least exp(A) — 1. This 
is certainly greater than the A given in Equation ( fl^) . and can be significantly larger for 
large A. 



Prakasa Rao and Sreehari [13] show that the property Rx = Var {X) characterizes the 
Poisson distribution (up to integer shifts) and Srivastava and Sreehari extend this 
to a characterization via weighted Poincare constants. 



Remark 3.4 Equation / fl^) in Example \3.3\ can also be proved using Poisson-Charlier 
polynomials ci!'\x), under the normalisation that Y^'^^Qllx{x)ci^\x)cm\x) = n\X"'6mn- 
It is well known (see for example JTS']) that these polynomials behave well with respect 
to A, with Ac^\x) = nc^'^i(a;). This approach exactly parallels the calculation using 
Hermite polynomials performed by Chernoff ^ to determine the Poincare constant of 
the normal distribution. 



Next, we show that our method reproduces the bound imphed by Table 2.2 of Klaasen 
[9] for binomial random variables. (Our method was specifically intended to generalise 
the result Klaasen gave in this case). 



Example 3.5 For Px{x) the mass function of a binomial B{n,p) random variable, 

in)f ^ _ xPxjx) _ p 

Px{x-l){n-x + l~ 1-p 



In other words. Equation fT^) holds with equality for D = p, so Lemma l2~< \ii) implies 
that the binomial B{n,p) has Poincare constant < np. 

Comparing this with the lower bound np{l—p) from Lemma Wl\ we can see this bound is 
reasonably accurate for a range of parameter values, certainly compared with the 1/(1 — 
p)" — 1 that follows from Bobkov and Gotze's Theorem \3.2[ 
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These calculations can presumably be performed in the spirit of Remark \3. 41 using the 
polynomials orthogonal with respect to binomial weights derived by Young in fW^). 



Klaasen's method, and the similar results of Cacoullos and Papathanasiou [S], in general 
give a weighted Poincare constant, while we prefer to consider the standard unweighted 
Poincare constant. 

Example 3.6 If X is the sum of n independent Bernoulli Bern (pj) random variables, 



then EX = Yll=iPi ^'^'^ (^) = Yl^=iPi{^ ~ Pi); so Theorem\Ll\i) gives Rx < 



Sr=iPj + l/2 + + Yll=iP'i (tighter bounds can be found using Theorem M.^ii) , but 



they are less transparent) . In this case, Klaasen does not give a bound on Rx, and the 
upper bounds implied by Bobkov and Gotze's Theorem \3.2\ are again not tight. 

Lemma \3. 1\ gives a lower bound Rx > ^"=i(Pi ~p'i)> so that upper and lower bounds 
are asymptotically within 1 of each other in the regime where Yl^pf tends to zero. It is 
perhaps not surprising that Le Cam flW proved Poisson approximation bounds in this 
regime. 

In fact, tighter upper bounds can be found, using the following result: 

Lemma 3.7 If X and Y are independent then Rx+y < Rx + Ry- 

Proof We consider g such that ^g{X + Y) = 0, and define h{u) = KYg{u + Y), which 
has the property that E,h{X) = 0. Now, as in the proof of Theorem 2(vii) of j^: 

Eg^{X + Y) = Ex [Eg'^iX + Y)\X] 

= Exye.T{g{X + Y)\X)+Ex[Eg{X + Y)\Xf 

< RyEx [Ey(A^)(X + Yf\X] + EhiX)"^ 

< RYE{Ag){X + Yf + RxE{Ah){Xf. 

For each x, 

Ah{x) = h{x + 1) - h{x) = J2 My) {9{x + y + l)-g{x + y)) , 

so using Cauchy-Schwarz: 

{Ah){xf < My) {9{x + y + l)-g{x + y)f , 

y 

and taking expectations using weights Px{x)^ we deduce that E(A/i)(X)^ < E{Ag){X + 
F)^, and the result follows. □ 
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Example 3.8 As in Example \3.5l take X to be the sum of independent Bernoulli 
Bern (pi) random variables. Since Example \3.^ shows that each Bern (pj) has i?Bem{pi) = 



Pi, using Lemma 3/l_ we deduce that Rx < X]r=iP« 



The improvement to the resuhing bounds is not surprising, since for example ^ shows 
that tensorization provides good bounds. However, we believe that the fact that The- 
orem 11.21 holds in general, for cases where such tensorization is not possible, makes it 
worthwhile despite the weaker bounds it produces in these special cases. [Note the dif- 
ferent value of the Poincare constant of Bern {pi) given in [2] is due to their restricting 
attention to functions g with g{2) = g{0)]. 

In fact, in the spirit of [8], it should be possible to give stronger results than Lemma 
13.71 in the case where X and Y are not Poisson, and hence improve the bounds given in 
Example 13. 8[ We do not pursue this route here. 



4 Remark 

We briefly mention a mistake in a previous paper: 

Remark 4.1 Part (ix) of Theorem 2 of /5/ considers mixtures of continuous random 
variables, defining by P(Z„ < z) = aP(X < z) + (I - a)¥iY < z). Although J^/ 
claims that the (continuous) Poincare constant satisfies Rz^ < max(i?X; Ry), this result 
is not true, in fact there is a sign error in the proof. 

The key observation is that if Rx < oo then the support of X must be an interval. 
Considering X and Y to be uniform on [0, 1] and [3, 4] respectively, J^/ shows that Rx 
and Ry are finite. However Rz = oo, since we can take functions arbitrarily close to 

, , f — (1 — a) for X < 1 
^ ' y a for x > 6 

which is non-constant on the support of Z^, but has zero derivative on the support of the 
mixture. 

A similar argument works in the N-valued case by taking X and Y uniform on {0, 1} 
and {3, 4} respectively. 
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